Thermal lensing-induced bifocusing of spatial solitons in Kerr-type optical media 
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Thermo-optical effects cause a bifocusing of incoming beams in optical media, due to the birefrin- 
gence created by a thermal lens that can resolve the incoming beams into two-component signals 
of different polarizations. We propose a non-perturbative theoretical description of the process of 
formation of double-pulse solitons in Kerr optical media with a thermally-induced birefringence, 
based on solving simultaneously the heat equation and the propagation equation for a beam in a 
one-dimensional medium with uniform heat flux load. By means of a non-isospectral Inverse Scat- 
tering Transform assuming an initial solution with a pulse shape, a one-soliton solution to the wave 
equation is obtained that represents a double-pulse beam which characteristic properties depend 
strongly on the profile of heat spatial distribution. 

PACS numbers: 42.65.Tg, 42.70.Gi, 44.10.-|-i 



In dielectric media with thermo-optical effects, the 
modulations of incoming beams can lead to a wavefront 
distortion reflecting their instability. Generally 

this instability gives rise to a depolarization of a high- 
power field d-Q due to a thermally-induced birefringence 
which is attributed [l0l-IT2| to a change in the refrac- 
tive index of the medium. It has been established [ij 
that this change in the refractive index originates from 
heat deposition in the propagation medium, resulting in 
a space-dependent temperature gradient (so-called ther- 
mal lensing) [l], 0, For media with linear indices, 
the thermal lens leads to a drastic change in the irradi- 
ance along the beam axis so that the resulting depolar- 
ization can strongly degrade the beam quality requiring 
thermal lensing compensation. However in nonlinear me- 
dia such as Kerr media, the nonlinearity can be a rele- 
vant self-compensaton factor [l^ stabilizing multi-wave 
modes generated by the thermal birefringence. 
While several materials exhibiting thermo-optical effects 
are known in the literature [l|, a most investigated one is 
the solid-state laser NdiYAG This material 
is represented as a rod crystal with a cylindrical geom- 
etry, where the change in temperature induces thermal 
distortion of incoming laser beams Q . For this particular 
material, several theoretical attempts have been made to 
formulate the spatial profile of the temperature gradient 
along the rod exploiting available experimental data. In 
particular, in refs. [l|, lisl [Tsj it was found that in the 
cylindrical rod configuration where the heat is generated 
at a constant rate a quadratic spatial distribu- 
tion gives a very good description of the experimentally 
observed birefringence and the resulting beam bifocus- 
ing 

But thermo-optical processes are actually common to a 
broad class of materials, not just solid-state lasers. In- 
deed, photonic crystals and optical fibers (including laser 
fibers) displaying nonlocal thermal and photothermal 
properties have been considered in the recent past, from 
both experimental and theoretical points of view [13, [3 ■ 
These materials share in common the fact that heat re- 



sulting from the input pump source causes physical vari- 
ation of the material. Namely, the material expands 
with the heat load due to a stress gradient formed which 
produces space-dependent birefringence in the material. 
Quite remarkably when the thermal gradient is strong 
enough the thermally-induced birefringence can resolve 
a polarized high-power input beam into two-component 
beams 0, S Il9| |- moreover upon recombination after 
traversing the bulk the two beams with two polarization 
directions are no longer in phase with one another such 
that the polarization state of the input beam cannot be 
recovered [3. 19.]. 

The thermally-induced depolarization phenomenon to- 
gether with the resulting beam bifocusing have been in- 
vestigated experimentally and theoretically by several 
authors ■ photonic crystals with a nonlinear Kerr- 
type dielectric susceptibility competing with diffraction, 
nonlocal thermal properties of the propagation medium 
have very recently been shown to induce twin-mode spa- 
tial solitons so-called dipole solitons (20l - |23j . The robust- 
ness of these double-pulse structures provides strong in- 
dication that the competition between nonlinearity and 
thermal lensing might be a stabilizing factor for the 
double-polarized modes, potentially observable in several 
optical materials including solid-state and optical fibers 
with thermally-induced birefringence. 
The present work aims at proposing a non-perturbative 
description of the generation of a double-polarized pulse 
beam in optical materials with thermally-induced lin- 
ear birefringence, paying particular attention to non- 
linear optical media with a quadratic spatial profile of 
temperature distribution. Our objective is to provide a 
consistent theoretical framework for understanding bet- 
ter the effects of thermal lensing on beam propagation 
in nonlinear optical materials, of which the YAG fiber 
laser which temperature profile was exprimentally estab- 
lished [U, m, We will first demonstrate, by solving 
the heat equation, that the topography of temperature 
distribution is quadratic in the case of a one-dimensional 
(ID) anisotropic material uniformly loaded at a constant 
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heat flux rate. Next, using the relationship between the 
temperature gradient and the induced inhomogeneous 
optical index, we shall derive an effective refractive in- 
dex for the material by combining the thermally-induced 
optical index and a Kerr-type refractive index reflecting 
the intrinsic nonlinearity of optical properties of the ID 
anisotropic material. With the help of this effective re- 
fractive index we shall formulate the propagation equa- 
tion for beams in the thermal nonlinear medium. Note 
that this equation is a Nonlinear Schrodinger equation 
(NLSE) with a repulsive external quad ratic potential, 
that has already been derived in ref. [2J1 but solved pcr- 
turbatively. Here instead, a non-perturbative treatment 
will be proposed based on a non-isospectral Inverse Scat- 
tering Transform (1ST) method [IM UHl with emphasis 
on IST's initial solution being single-pulse shaped. 
Thermal lensing can be described in simple words as a 
thermo-optical process associated with a weak absorp- 
tion of an input beam that induces a nonzero temperature 
gradient across the material, leading to a spatial variation 
of its refractive index [l|, [l^ [2l| . Recent experimental as 
well as theoretical developments on this process suggest 
that the underlying mechanism involves a local refractive 
index change An{x) which increases linearly [H,[2l| with 
the temperature change AT{x) i.e. An — (3AT, where 
/? = dn/dT refers to the thermo-optic coefficient |2TI.[27|. 
Thus when the optical beam of a uniform thermal load 
gets slightly absorbed and heats the material, this pro- 
duces heat that is conveyed by the electromagnetic wave. 
If p denotes the uniform heat flux density and k the heat 
conductivity coefficient, the heat diffusion in the material 
along a preferred direction (for ID anisotropic materials 
of current interest), driven by the uniform heat flux load, 
is determined by the heat equation: 

KV'^T{x)^-p. (1) 

Since eq. ([IJ is key to the current analysis we consider 
its most general solution given by: 

T{x) ^ -^a2x'^ + aix + To, (2) 

where 02 — p/n, ai and Tq being two arbitrary real con- 
stants. Formula ^ is consistent with the quadratic law 
of temperature variation found for most laser fibers in 
the presence of a uniform head load. More specifically 
in YAG fiber lasers this law is common [l|, [l^, [l^ and is 
consistent with the optical bifocalility associated with a 
thermally-induced birefringence, that promotes double- 
polarized laser beams from an input laser field. 
For the sake of simplicity we require the temperature gra- 
dient to be zero and temperature to take a bare value Tq 
at X = (i.e. the ambient temperature). The change of 
temperature in the material along the x axis, hereafter 
assumed to be the axis of beam propagation, then reads: 

AT(x) - To - Tix) 

= -jaix^. (3) 



With the last formula we derive the following expression 
for the local refractive index change: 

An{x) = — aa;^, a = pP/n. (4) 

Now if the intrinsic optical properties of the material 
are dominated by Kerr-type phenomena, the homoge- 
neous part of the refractive index can be expressed as 
(fc2c/27r)n2 / where / is the beam intensity. With the 
help of (HI we can readily define an effective refractive 
index n{x) for the thermal nonlinear material viz: 

n[x, I] An{x) + {k^c/2TT)n2 I. (5) 

Assuming that the wave motion is fast along the axis 
of anisotropy (i.e. x) but very slow along z [2l|, the 
paraxial approximation on the 2D wave equation for an 
electromagnetic field q{x, z) leads to: 

g+2zfcg + ,[,,/],^0. (6) 

As already stressed eq. dH) is actually not new, in- 
deed the same equation was obtained [2J] for the same 
problem but solved following the collective-coordinate 
method. To this last point, eq. (jG]) is an inhomoge- 
neous NLSE and so can in principle be solved using the 
collective-coordinate method. However this is a pertur- 
bative method and consequently requires that the ther- 
mal lensing is sufficiently weak, so that the Kerr nonlin- 
earity remains the main governing factor in the modu- 
lation and stability of signals in the thermal nonlinear 
medium. So to say any input beam sent in the medium 
must be modulated into a signal of permanent single- 
pulse shape, with eventually an acceleration or slow down 
of the pulse due to the thermal lensing. In fact this 
consideration is very far from any acceptable consistency 
with the physics of the process under study, which the 
double-polarization of the incident beam is a most salient 
aspect. 

Being interested in a solution to eq. ([6]) which is more con- 
sistent with experiments we try for a non-perturbative 
approach. Remark to start that this equation can be 
rewritten in the following form: 

+ 2»fc|^ + [{ec/2n)n2 \q\'' - V{x)] q = 0, (7) 

corresponding to a NLSE with an external potential 

V{x)^~ax^/2 (8) 

which is quadratic in x with a maximum at x = 0. The 
physics behind this quadratic potential is contained both 
in its profile and the parameter a defined in Q , which de- 
termines the strength of the thermal birefringence on the 
beam shape. One remarkable side of this physics emerges 
from the assumption of the heat fiux density p and the 
heat conductivity k as being fixed, such that a appears to 
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be increasing with an increase of the thermo-optic coef- 
ficient /3. Thus when /? is increased the curvature of the 
quadratic profile of heat distribution in the material is 
more and more pronounced so that the effect of thermal 
lensing on beam modulation is more and more strong. In 
fact the external potential is expected to be much effec- 
tive on the beam position and according to the profile 
of this potential given by (|5]), the centre of mass of the 
input beam should experience a trapping force from the 
manifestly expulsive quadratic potential. 
When a is large enough such that the potential field 
erected by the thermal lensing process on the beam path 
is strongly localized, its contribution must be fully taken 
into consideration. In this last respect we follow the non- 
isospectral 1ST technique proper '25', '26'] to equations of 
this specific kind, considering an initial signal q{x, z = 0) 
of a permanent single-pulse shape along a; at z = 0. For 
eq. ^ this technique leads to the following one-soliton 
solution: 



(a) 



(a) 



with 



q{x, z) = Qo{x, z) sech [^^(x, z)] e**(^'^\ 



Qo{x,z) = J 

V cn2 k 



^{x, z) 
$(a;, z) 



ln(2|/(z)|) 
2V{x)g{z) - 



- 2V{x)fiz) 

$0, 



sech?{z/Q 
-8(p/3/k)2 tanh2(0/C)' 



tanh(z/C) 



2kV2^^ 1-|-8(p/3/k)2 tanh2(z/C)' 
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If the "sech" function in ([9]) reminds a pulse signal, the 
complicated dependence of its prefactor Qo{x^ z) in x and 
z clearly suggests an actually complex pulse structure for 
the beam on the xz plane. To gain insight about what 
this dependence implies for the signal profile, in fig. [T] we 
plotted the beam amplitude |g(a;, z)! as a function of x 
at two different positions in the direction z transverse to 
the beam propagation, and for four distinct values of a. 
More explicitely the four left graphs in fig. [T] represent 
|g(x, z)! versus x z = 0.5 for a = 0.1, 0.05, 0.01 and 
0.005, while the left graphs represent [(/(x, z)| versus x 
for same values of a but at z = 10. As one sees, the sig- 
nal intensity is a two-component pulse which intensities 
are strongly dependent on the strength of the thermal 
lens potential. It is quite noticeable on exploring the 
six figures that when a increases, the intensities of the 
two-component pulse increase while their width at half 
tails diminish. The last behaviour is consistent with the 
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FIG. 1: a:-£ixis profiles of \q{x, z) \ for z = 0.5 (left graplis) and 
z = 10 (right graphs). Values of a are (a) 0.1, (b) 0.05, (c) 
0.01, (d) 0.005 



a dependence of parameter C, defined in (1151) . which in- 
deed represents the average spatial extension of the pulse 
along the x axis. 

Another relevant feature emerging from graphs of fig. [D 
is the fact that when a is decreased for a fixed value 
of z, the two constituents pulses in the double-polarized 
beam preserve their shapes but their peak positions (i.e. 
centres) are gradually shifted. The last feature is more 
transparent in the left graphs corresponding to a rela- 
tively large value of z. In fact, the last behaviour can be 
interpreted in terms of a ring profile for the signal inten- 
sity in the xz plane as refiected by the contour plots of 
fig. [5J where shadows of the double-pulse soliton in the 
xz plane are represented for different values of a. 

The figure clearly indicates an increase of the separa- 
tion between pulses in the double-pulse signal, implying 
that the radius and curvature of the ring signal in the 
xz plane are fixed by the magnitude of a. 
The double-pulse structure obtained as well as its 
ring profile emerging in 2D are reminiscent of dipole 
and ring-vortex solitons, predicted recently in some 
nonlocal nonlinear media [2l|, [l^, [2^ . However in these 
previous studies the physical origins of nonlocalities 
were generally not well specified, whereas in our context 
the quadratic form of the optical inhomogeneity has an 
experimental foundation [ll[ll,[ll]. Also, a previous 
attempt [l^l to model the same problem led to eq.®- 
In this previous work the equation was treated pertur- 
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batively and so results could not reflect the remarkable 
aspects underlying the physics of thermal lensing in 
Kerr media. As end remark, the universality of the 
1ST one-soliton solution ^ can be checked by applying 
any other exact spectral method to the generating 
equation, such as the Darboux method 29] combined 
with a Lax-pair formalism with non-conserved spectral 
parameters [30ll3l|. 
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FIG. 2: Contour plots of |q(a::,z)| for different a listed as: (a) 
0.01, (b) 0.005, (c)O.OOl, (d) 0.0005, (e) 0.0001, (f) 0.00001. 
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